A Precision Measurement of the Muon Decay Parameters p and 5 
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The TWIST collaboration has performed new measurements of two of the parameters that describe 
muon decay: p, which governs the shape of the overall momentum spectrum, and S, which governs 
the momentum dependence of the parity-violating decay asymmetry. This analysis gives the results 
p = 0.75014 ± O.Q0017(stat) ± Q.00044(syst) ± 0.00011(r;), where the last uncertainty arises from the 
correlation between p and the decay parameter rj, and 5 = 0.75067 ± 0.00030(stat) ± 0.00067(syst). 
These are consistent with the value of 3/4 given for both parameters in the Standard Model of 
particle physics, and are a factor of two more precise than the measurements previously published 
by TWIST. A new global analysis of all available muon decay data incorporating these results is 
presented. Improved lower and upper limits on the decay parameter P^^ of 0.99524 < P^^ < 
^ < 1.00091 at 90% confidence are determined, where P^ is the polarization of the muon when 
it is created during pion decay, and ^ governs the muon decay asymmetry. These results set new 
model-independent constraints on the possible weak interactions of right-handed particles. Specific 
implications for left-right symmetric models are discussed. 

PACS numbers: 13.35.Bv, 14.60.Ef, 12.60.Cn 



I. INTRODUCTION 

In the Standard Model (SM) of particle physics, 
the charged-current weak interaction violates parity 
maximally — only left-handed particles (or right-handed 
antiparticles) are affected. The TWIST experiment is 
a high-precision search for contributions from non-SM 
forms of the charged-current weak interaction, including 
parity-conserving currents. 

The decay of the positive muon into a positron and two 
neutrinos, /i^ — » e^i>i,v^, is a purely leptonic process, 
making it an excellent system for high-precision studies 
of the weak interaction. It proceeds through the charged 
weak current — mediated by the W boson — and can be 
described to a good approximation as a four-fermion 
point interaction. The matrix element for the most gen- 
eral Lorentz-invariant, local, four-fermion description of 



muon decay can then be written as 
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where the gju specify the scalar, vector, and ten- 
sor couplings between /i-handed muons and e-handed 
positrons [1], and satisfy certain normalizations and con- 
straints. In the Standard Model, g^j^ — 1 and all 
other coupling constants are zero. The probability Q^^ 
(e, ^ = L, R) for the decay of a ^-handed muon into an 
e-handed positron is given by 
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where 6^^, — 1 for e — ^ and Stfi — for e ^ [i. The 
probability: 
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sets a model independent limit on any muon right-handed 
couphngs P, Q- 



The differential muon decay spectrum [3, |j, Q can be 
described following the notation of Fetscher and Ger- 
ber j^] as 



d^r 



dxd{cos0s) Air^ 



■A' 7714 ^2 






(4) 



X {Fis{x)+Pf,^ COS 9,Fas{x)) 



where Gp is the Fermi coupling constant, 9s is the an- 
gle between the muon spin and the positron momentum, 
Emax ~ 52.8 MeV is the kinematic maximum positron 
energy, x = Ee/Emax is the positron's reduced energy, 
Xq = rUe/ Emax is the minimum possible value of x, cor- 
responding to a positron at rest, and P^ is the degree 
of muon polarization at the time of decay. P^ can be 
used to determine P^, the degree of muon polarization 
at its creation from pion decay, when the amount of de- 
polarization undergone by the muon is known. The two 
components of the decay spectrum in Eqn. Q are the 
isotropic component: 
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and the anisotropic component: 
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R.C. represents the electromagnetic radiative corrections, 
which have been calculated to 0{a^). Corrections due to 
the strong interaction in loops give a fractional contri- 
bution on the order of 4 x 10^^ [7], which is more than 
two orders of magnitude smaller than the ultimate pre- 
cision goals of TWIST. The quantities p, S, r], and ^, 
often called the Michel parameters, are bilinear combi- 
nations of the weak coupling constants and describe the 
shape of the decay spectrum. These decay parameters 
can be used in combination with other muon decay mea- 



surements, such as inverse muon decay (e 



A* f^e) 



and the polarization of the decay positron, to determine 
limits on the weak coupling constants. 

Left-Right Symmetric (LRS) models [Sj] comprise an 
interesting class of extensions to the SM. These mod- 
els include a right-handed weak coupling, which is sup- 
pressed by the mass of the associated gauge boson. LRS 
models contain four charged gauge bosons (Wi , W2 ) 
with masses mi and 7712, and two additional massive neu- 
tral gauge bosons. The mass eigenstates Wi and W2 are 
related to the weak eigenstates Wl and Wr through a 
mixing angle C.. ql and gn are the coupling strengths 
of the LRS weak interaction to left- and right-handed 
particles. LRS models affect the muon decay spectrum. 



including a modification to the decay parameter p: 
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therefore requiring p < 0.75. To a good approximation, 
the value of S is unaffected by LRS models. 

Many other proposed SM extensions also lead to mod- 
ifications of the Michel parameters. For example, super- 
symmetric models can lead to a non-zero value for 5^^ 
|9|. The values of ^ and 6 can be combined to provide 
a model-independent limit on right-handed couplings in 
muon decay: 
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Previous measurements by TWIST had set new lim- 
its, oi p = 0.75080 ± 0.00032(stat.) ± 0.00097(syst.) ± 
0.00023(7?) [13] and (5 = 0.74964 ± 0.00066(stat.) ± 
0.00112(syst.) [Il|. This paper presents new results 
from a refined analysis of newer data, providing details 
about the TWIST experiment and analysis (Sects. HIl 
and IIV[) including improvements over the previous stud- 
ies (Sect.|V|, discusses the validation of the Monte Carlo 
simulations (Sect. lyi]) . describes the current state of the 
systematic uncertainties in the experiment (Sect. IVIip . 
and presents new measurements of p and S (Sect. IVIIip . 
These new results are then incorporated into a global 
analysis of all muon decay data (Sect. IIX|) . 



II. EXPERIMENTAL PROCEDURES 
A. Experimental Setup 

Highly polarized muons resulted from the decay of pi- 
ous stopping at the surface of a carbon production target 
bombarded by 500 MeV protons at TRIUMF. The M13 
beamline |12l | selected positive muons with a momentum 
of 29.6 MeV/c, with a momentum bite of Ap/p w 1% 
FWHM, and delivered these in vacuum to the TWIST 
spectrometer [lj| at a typical rate of 2.5x10'^ per sec- 
ond. The beam had a contamination of positrons, at the 
same momentum, with a typical rate of 22 kHz, as well 
as a small fraction of pions. 

The muon beam was characterized and tuned using a 
low pressure (8 kPa dimethylether (DME) gas) remov- 
able beam monitoring chamber located upstream of the 
TWIST spectrometer 



141. The beam monitor was in- 



serted for measurement of the beam properties, and re- 
moved during data-taking. 

The TWIST spectrometer fl3l was designed to mea- 
sure a broad range of the muon decay spectrum, allowing 
the simultaneous determination of the decay parameters. 
The detector consisted of 44 drift chambers (DCs) and 
12 multiwire proportional chambers (PCs) in a planar 
geometry, symmetric about a muon stopping target foil 



at the centre. Figure [T] shows the upstream half of the 
detector and the four PCs surrounding the stopping tar- 
get. The chambers were placed in a highly uniform 2 T 
solenoidal magnetic field. The z axis was defined to be 
the detector axis. 
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FIG. 1: (Color online.) Side view of the upstream half of 
the TWIST spectrometer planar chambers and support struc- 
tures. Muons stopped in the target foil, which also served as 
a chamber cathode. The use of precision spacers and tension 
control ensured that the chamber wire positions were known 
to 30 /im in z. The spectrometer is symmetric about the 
target foil. 

To reduce scattering and to allow the muons to reach 
the central target foil, the detector was designed to be 
very thin; there was approximately 140 mg/cm^ of ma- 
terial from the vacuum of the M13 beamline through to 
the centre of the stopping target. 

The muon stopping target was a 71 ± 1 /im foil of 
99.999% pure aluminum, which also served as a cathode 
foil for the adjacent two PCs (see below). The muon 
stopping distribution, as determined from the last DC 
or PC plane in which the muon left a signal, was used 
in a feedback loop to control the fractions of He and 
CO2 in a gas degrader at the end of the beamline, in 
order to maintain the average stopping position of the 
selected muons at the centre of the target. Muons were 
required to be recorded by the PC immediately before the 
stopping target (PC 6) and not by the PC immediately 
after the target (PC 7); simulations showed that 97.0 ± 
1.5% of selected muons stopped in the target, with the 
rest stopping in the CF4/isobutane, the cathode foils, or 
the wires in the vicinity of the target. The experiment 
watched for muons decaying at rest into positrons; decays 
of muons in flight were identified and discarded. 

Each PC and DC consisted of a wire plane and two 
cathode foil planes, all oriented perpendicular to the de- 
tector axis. Chambers were grouped together in modules 



as described below, and the cathode foils interior to each 
module were shared between adjacent chambers. Each 
wire plane was in either U ot V orientation, and these 
were at right angles to each other and 45° to the vertical. 

Four PCs were located at each end of the detector, 
and two on either side of the stopping target. Each PC 
included 160 sense wires at 2 mm pitch. A fast drift gas 
(CF4/isobutane) was used in these chambers. The PCs 
mainly provided particle timing information; in addition 
the widths of the timing signals from the PCs were used 
to discriminate muons from positrons. 

Each half of the detector included 22 DCs, with a slow 
drift gas (DME), which had a small Lorentz angle; these 
provided precise measurements of the e "*" p osition as it 
crossed the chamber. Each DC plane [1^ included 80 
sense wires at 4 mm pitch. The DC planes were very 
slightly asymmetric in the direction along the detector 
axis, in that the wires were located 150 iim. off of the 
centre of the cell. 28 of the DCs were paired together in 
14 modules, with U and V wire planes in each pair. Ad- 
ditional "dense" stacks of eight planes each were located 
near either end of the detector. 

The space between the chambers was filled with a 
3% N2 and 97% He gas, minimizing the material thick- 
ness of the detector. The chamber gas and the helium 
mixture were maintained at atmospheric pressure. The 
differential pressure between the chambers and the sur- 
rounding volumes was controlled to stabilize the positions 
of the chamber foils. 

Chamber alignments perpendicular to the detector axis 
were measured using 120 MeV/c pion tracks with the 
magnetic field off. Alignment measurements were taken 
at the beginning and end of the data-taking period. 
These data were also used to measure wire time offsets 
introduced by the electronics and cabling. The plane 
positions and rotations about the beam direction were 
determined to an accuracy of 10 /im and 0.4 mrad. Rel- 
ative wire positions were known to 3 /^m. The alignment 
of the chambers to the magnetic field was measured using 
positron tracks in the 2 T field, and had an uncertainty 
of 0.03 mrad. 

The chambers were positioned within the 1 m bore 
of a liquid-helium-cooled superconducting solenoid. The 
magnet was placed inside a cube-shaped steel yoke, ap- 
proximately 3 m per side, designed to increase the unifor- 
mity of the field within the detector region. The shape 
of the z component of the magnetic field was mapped 
using an array of Hall probes mounted on a radial arm, 
aligned along the solenoid axis and calibrated using an 
NMR probe. The probes' positions at each field sample 
were known to ±1 mm, and the field was mapped to a 
precision of ±1 x 10^"* T. The standard operating setting 
for TWIST is 2 T at the centre of the solenoid; at this 
setting the field is uniform to within 8x 10"^ T within the 
tracking region. A finite-element simulation of the mag- 
netic field was performed using Opera3d "It?!, provid- 
ing the full three-dimensional magnetic field throughout 
the detector volume and well outside the measured re- 



gion. Within the tracking region the simulated field map 
agreed with the measurements to within ±2 x 10~^ T. 

Event acquisition was triggered by a thin scintilla- 
tor upstream of the spectrometer. Only 1.8% of events 
were triggered by a beam positron, in spite of the much 
higher rate of beam positrons compared to muons; most 
beam positrons were discriminated against by the trig- 
ger threshold. Remaining beam positrons were identified 
with high efficiency in the data since they left signals in 
the full length of the detector. 

Preamplifier chips were mounted directly on the cham- 
bers and drove custom postamplifier and discriminator 
modules placed 2 m away. The chamber behavior was 
exceptionally stable, with less than one high voltage trip 
per month. All output channels were functional. 

Data acquisition from scintillators and tracking cham- 
bers [17| was performed using LeCroy Model 1877 TDCs. 
The trigger and TDC read-out recorded signals in 0.5 ns 
time bins from 6 fis before to 10 fis after a muon passed 
through the trigger scintillator. The start and stop of 
each pulse was recorded so that the pulse width could be 
determined; these widths can be related to the amount 
of energy deposited in the cell. In this configuration up 
to eight wire signals could be recorded for each wire in 
any triggered event. A fixed blanking time of 80 /iS was 
imposed after each accepted trigger, to allow each TDC 
to finish conversion before the next event was recorded. 

The data considered for this analysis were more than 
1.5x10^ muon decays, taken during 2004. The data sets, 
summarized in Table|Tl were taken under a variety of con- 
ditions (low polarization from beam steering, rate, muon 
stopping position, etc.). Provided the simulation repro- 
duced these conditions correctly, the decay parameters 
extracted from the data should be independent of the 
run conditions. 

TABLE I: Description of data sets, in dironological order. 



Description 



Events Accepted 
(xlO^) (xlO^) 



Low Polarization, 


Centred 


209 


8 


Low Polarization, 


PCS Stops 


94 


2 


Centred Stops 




287 


11 


3/4 Stops A 




323 


12 


High Rate 




198 


7 


Aperture 




263 


9 


3/4 Stops B 




157 


6 



Total 



1531 



55 



III. SIMULATION 

The decay parameters were extracted from the data by 
comparing the data to a simulated spectrum, as discussed 
in Sect. lIVDl Monte Carlo (MC) simulations were run to 
match the conditions of each of the seven main data sets. 



in addition to simulations run for studying systematic un- 
certainties. Each simulation included 2-3 times as many 
muon decays as the corresponding data set. The simu- 
lations were implemented using Geant 3.21. The output 
was in exactly the same format as the files produced by 
the data acquisition system, and were processed in the 
same way as real data. Space-time relationships (STR) 
determined with Garfield [1^ were used to model drift 
chamber response. 

The theoretical decay spectrum included full radia- 
tive corrections at 0{a) [T^, as well as O(a^-L^) and 
0{a^L) [13, [m, where L = \og{ml/ml) « 10.66. 
0{a^L°) terms have also been calculated [22; the ef- 
fect of neglecting these last terms has been evaluated and 
shown to be negligible (see Sect. lVIIFJ) . All radiative cor- 
rections are calculated within the Standard Model. The 
values of the muon decay parameters assumed by the 
simulation in its theoretical decay spectrum were kept 
hidden until the end of the study; see Sect. IIVDI below 
for details. 

The simulation included all known depolarization ef- 
fects, including interactions with the magnetic field as 
the muon enters the solenoid and depolarization in the 
stopping target. 



IV. DATA ANALYSIS 

The full analysis procedure was applied in the same 
way to both data and simulation. To the level that the 
simulation accurately represented the data, this canceled 
spectrum distortions due to detector response, positron 
energy loss and scattering, reconstruction biases, and 
other effects, which would otherwise lead to systematic 
errors in the measurement. 

Due to the large amount of simulation and analysis 
required, the Western Canada Research Grid (WestGrid) 
was used. TWIST used approximately 10,000 CPU-days 
for this simulation and analysis. 



A. Event Reconstruction 

To reconstruct an event, the hits — signal times on in- 
dividual wires — were first grouped based on timing infor- 
mation from the PCs. Tracks within these groups were 
then identified using the distribution of DC hits in space 
and time, as well as the hit widths. DC hits associated 
with decay positrons were then used to reconstruct the 
energy and angle of each positron. 

Pattern recognition was performed on the decay 
positrons using the spatial hit distributions, to determine 
an initial estimate of the positron track. The track fit pa- 
rameters were the position and momentum three- vectors 
at the DC closest to the target; the main parameters of 
interest were the positron's total momentum p, and the 
angle 9 between the positron momentum vector and the 



detector axis. Multiple overlapping tracks could be dis- 
tinguished at this stage. The initial track estimate was 
then refined in a x^ fit using the hits' drift times, in com- 
bination with maps of the STRs of the drift chamber cell 
as provided by the Garfield chamber simulation soft- 
ware 18]. The positron track was assumed to deviate 
from a helix by continuous energy loss through the gas 
volumes and discrete energy loss through each foil en- 
countered, using mean energy loss formulas [S]- Tracks 
were allowed discrete deflections at each DC module and 
in the dense stacks. The deflection angles were flt pa- 
rameters, with associated penalties to the fit x^, based 
on the method described by Lutz [23| . 

Event and track selection cuts were then applied and 
the decay spectrum was assembled. A muon hit was re- 
quired in PC6, immediately upstream of the stopping tar- 
get, to ensure the muon reached the target. A hit in PC7, 
immediately downstream, vetoed the event. Events with 
multiple muons were rejected. The muon was required 
to decay between 1 /zs and 9 fis after the trigger; the de- 
lay ensured that the ionization from the incident muon 
was collected before a hit from the positron was recorded. 
Events where a beam positron arrived within 1 /is of ei- 
ther the muon or decay were rejected as well. Additional 
cuts included the muon flight time through the M13 beam 
line, used to reject muons from pion decays in flight, and 
a requirement that the muon stopped within 2.5 cm of 
the detector axis, which ensured that all decay positron 
tracks within the fiducial region (Sect. IIVC[) were fully 
contained within the detector. Note that these cuts de- 
pend on observations of the muon prior to decay or of 
the time of decay, and are hence unbiased in terms of the 
properties of the positron. 

Figure [5] shows a muon decay spectrum reconstructed 
from the "Centred Stops" data set (3x10^ muon decays; 
see Table HI). 
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FIG. 2: (Color online.) Reconstructed muon decay spectrum. 
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B. Energy Calibration 

Differences in the details of the muon stopping posi- 
tion within the target foil, differences in the STRs, the 
accuracy of the magnetic field maps, and knowledge of 
the wire positions resulted in differences in the energy 
scales of the simulation and the data. The kinematic 
endpoint region (near 52.3 < p < 53.4 MeV/c) was used 
to determine the relative energy calibration between the 
data and the simulation after reconstruction was com- 
plete. The relative positions of the endpoint were com- 
pared between data and simulation in bins of angle (see 
Fig. [3] for an example), and the differences were parame- 
terized separately for upstream and downstream parts of 
the spectrum, according to 



Ape = B, 



_Ai_ 
cos^' 



i = US, DS, 



(9) 



where Ape is the difference in the momentum of the spec- 
trum endpoint between data and simulation. Equation 



FIG. 3: (Color online.) Uncalibrated data and simulation 
histograms in the momentum range of the decay spectrum 
used for energy calibration, for —0.53 < cos6 < —0.50. 



[9] describes an energy loss proportional to the amount of 
material encountered by a particle; for a planar detector 
the amount of material is proportional to l/cos0. This 
form was found to describe Ape well, also. Uncertainties 
in the calibration parameters Bi and Ai were statistical. 
This relative calibration technique replaced the calibra- 
tion using an analytic approximation to the end point 
shape used in the previous analyses [10|, llH, 1241 • 



C. Fiducial Region 

Limitations of the reconstruction and physical aspects 
of the detector required the imposition of fiducial cuts, 
which were applied after the data and simulation were 
calibrated and assembled into spectra. Improvements 



in reconstruction allowed us to use a larger fiducial re- 
gion for this analysis than was used for previous TWIST 
measurements [Ifl, [ill, UM ■ The fiducial region adopted 
for this analysis required 0.50 < \cos6\ < 0.92, p < 
51.5 MeV/c, 10.0 < pt < 39.7 MeV/c, and \p^\ > 
13.7 MeV/c. This region is shown in Fig. |4l 
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FIG. 4: (Color online.) A reconstructed muon decay spec- 
trum; the dot density is proportional to the number of 
positron tracks reconstructed in a {p,cos6) bin. The outlines 
show the fiducial region used for this analysis. 



The 51.5 MeV/c momentum cut removed the region 
of the spectrum used by the energy calibration; this also 
removed the only sharp feature of the spectrum, greatly 
reducing the sensitivity of the decay parameters to the 
resolution. The minimum |cos0| cut elinrinated high an- 
gle events, which underwent large amounts of multiple 
scattering and were difficult to reconstruct. The maxi- 
mum I cos 1 cut was limited by increased reconstruction 
bias and a difference in reconstruction efficiency between 
data and simulation at small angles. The minimum pt 
cut eliminated tracks with radii too small to provide suf- 
ficient transverse separation in the hits for reliable fits. 
The maximum pt cut, in combination with the require- 
ment that muon stopped within 2.5 cm of the detector 
axis, ensured that all accepted decay positrons were fully 
contained in the detector. The minimum \pz\ cut re- 
moved tracks with wavelength comparable to a detec- 
tor periodicity of 12.4 cm. The robustness of the anal- 
ysis against variations in these cuts was checked using 
variations of 0.02 in cos 6* or 0.5 MeV/c for the momen- 
tum cuts; neither the decay parameter measurements nor 
the systematic uncertainties were changed by more than 
1 X 10-5. 

The fiducial cuts were applied to the spectrum his- 
togram, and were therefore influenced by the histogram 
binning. A histogram bin was included in the fit if the 
{p,cos9) value of its centre passed the above fiducial 
cuts. The widths of the spectrum histogram bins were 
0.5 MeV/c in momentum, and 0.02 in cos^. 



D. Spectrum Fitting 

The decay parameters were determined by compar- 
ing the shape of the two-dimensional data spectrum to 
that of a simulated spectrum with matching conditions. 
The muon decay spectrum is linear in the parameters 
{p,r],Pf^S,,Pfj,^S), so the data spectrum can be described 
in terms of the MC spectrum as 
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where Sd — d-T d / dx d{cos 6) is the decay spectrum 
from data, and Sm — d^T m / dx d{cos 6) is the spectrum 
simulated using Eqn. ([4]) to generate the muon decays. 
Note that the simulated spectrum is generated using hid- 
den values of the decay parameters. dS/dp etc. are the 
derivatives of Eqn. @ with respect to the decay param- 
eters. Radiative corrections were left out of these deriva- 
tives, under the assumption that the dependence of the 
radiative corrections on the decay parameters is negli- 
gible; the exception was dS/dP^S,, which included the 
anisotropic radiative corrections to facilitate the consis- 
tent treatment of P^^ as a product. 

"Derivative spectra" were fully simulated and analyzed 
spectra in the same way as Sm- To generate these spec- 
tra, the magnitudes of derivatives of Eqn. ([4]) with re- 
spect to the decay parameters {\dS/dp\ etc.) were treated 
as probability distributions for the purposes of spectrum 
generation; a sign was associated with each event accord- 
ing to the sign of the derivative at that point, and this 
sign was used as a weight when the reconstructed event 
was included in the final derivative spectrum. 

The coefficients Ap, A(P^^), and A(P^^(5) are the fit 
parameters, and represent the differences in the decay pa- 
rameters between data and simulation. The muon decay 
spectrum in our fiducial does not accurately determine 
T] because of the xq coefficient in Eqn. ([5]), so the value 
7] = —0.0036 ±0.0069 from the recent global analysis |25l | 
was assumed. 

Only differences between the real decay parameters in 
data and those assumed by the simulation were mea- 
sured directly in this procedure. The values of the 
parameters assumed by the simulation {ph,Sh,^h) were 
chosen randomly within 0.01 of the SM values while 
being constrained to physically meaningful values (e.g. 
£,S/p < 1), and these values were kept hidden throughout 
the analysis and the evaluation of systematic uncertain- 
ties. (5, and hence AS, was extracted from A(P^^(5) as 
i^hSh + A{P^^5))/{(h + A(i^^C))- When measuring sys- 
tematic uncertainties, it was sufficient to assume Stan- 
dard Model values for the extraction of AS. In this way 
the measurement outcome remained unknown until the 
values of the hidden parameters were revealed. 

The ability of the fitted Monte Carlo spectra to de- 
scribe the data can be viewed in terms two distributions: 



the angle-integrated momentum distribution of the decay 
positrons, which is governed by p (with 77), and the mo- 
mentum dependence of the decay asymmetry, governed 
by S. The momentum spectrum of accepted events from 
a standard data set is compared with that from the corre- 
sponding fitted simulation in Fig. [5] The decay asymme- 
try as a function of momentum for accepted events from 
the same data set is compared with the fitted simulation 
in Fig.[ni Here asymmetry is defined as {B — F)/{B + F), 
where B and F are the number of counts in the backward 
(cos 9 < 0) and forward directions in a given momentum 
bin. The normalized residuals demonstrate the quality 
of the fit and the lack of bias. 
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FIG. 5: (Color online.) Momentum distribution of accepted 
events, "Centred Stops" data set. Deviation from the shape of 
Eqn. Q below 28 MeV/c and above 45 MeV/c comes from the 
acceptance constraint shown in Fig. |4] The lower plot shows 
the difference between data and fit for each bin, normalized 
by the statistical uncertainty in the difference. 
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surement used more advanced analysis and simulation 
techniques than were available for the previous study, as 
described below. A new, independent set of hidden decay 
parameters was generated for use in the simulation, and 
this new set was used exclusively for this analysis. There- 
fore, as with the previous measurements from TWIST, 
the measurement outcome could not be known until the 
hidden decay parameters were revealed. 



V. SUBSTANTIAL CHANGES SINCE 
PREVIOUS STUDIES 

There were a number of improvements to the exper- 
imental apparatus and data-taking techniques between 
2002, when data were taken for the first TWIST mea- 
surements of p and i5 [l3,[ll|j and 2004, when the data 
were taken for both the first TWIST measurement of 
-P/T? 1231 and the present measurement of p and 5. The 
geometry and material of the muon stopping target used 
in 2004 was much better known than the target used 
previously. Moreover, additional monitoring and feed- 
back controls were implemented for the taking of data in 
2004. These improved the stability of the dipole mag- 
nets in the beamline, the fiatness of the chamber foils, 
average muon stop position, and other aspects of the ex- 
periment, significantly increasing the quality of the data. 
The new muon beam monitoring chamber allowed sig- 
nificantly better characterization of the muon beam for 
input into the simulation, as well. 

Since the previous measurements, the simulation and 
analysis were modified to account for the asymmetric 
construction of the drift chambers, where the wires were 
not centred between the foils (see Sect. Ill A1) . The track 
reconstruction used by the present analysis accounted for 
energy loss by the positrons, reducing reconstruction bias 
by several keV/c; other reconstruction improvements fur- 
ther reduced the reconstruction biases, particularly at 
low angles. Energy calibration used by the present anal- 
ysis was performed as a relative calibration between simu- 
lation and data, resulting in a more accurate calibration. 

As this analysis was nearing completion, a technique 
was developed to extract STRs directly from the data. 
This occurred too late to be incorporated into this en- 
tire analysis, but it was used to estimate the systematic 
uncertainty due to mismatch in the chamber response 
between data and MC (see Sect. I VII Bj) . 



FIG. 6: (Color online.) Decay asymmetry of accepted events, 
"Centred Stops" data set. The lower plot shows the difference 
between data and fit for each bin, normalized by the statistical 
uncertainty in the difference. 

This study used the same data analyzed by Jamieson 
et al. for the direct measurement of the decay parameter 
P^C, published in 2005 [13]. However, the P^£, analysis 
did not include the study of the p or 8 systematic uncer- 
tainties and corrections. Furthermore, the present mea- 



VI. SIMULATION VALIDATION 

Since the data are fitted with simulated spectra to mea- 
sure the muon decay parameters, the simulation must 
correctly reproduce physical and detector effects in an 
unbiased way; of particular importance are energy loss, 
multiple scattering, and the probabilities for the pro- 
duction of delta rays and bremsstrahlung. Direct com- 
parisons between simulation and data were performed 



with specialized data not used for decay parameter mea- 
surements. In this mode, labeled "upstream stops" , the 
momentum of the muon beam was reduced so that the 
muons stopped upstream of the DCs. Downstream-going 
decays then passed through the entire detector, and were 
reconstructed twice, first using only the upstream half of 
the detector, then using only the downstream half. En- 
ergy loss, scattering, helix fitter biases, and reconstruc- 
tion resolution all resulted in differences in the proper- 
ties of the two tracks. Distributions of Ap = pjjs — Pus, 
A6 = 9ds — Ous, and other differences were used to ex- 
amine reconstruction and physical effects, independent 
of the shape of the decay spectrum, and to compare the 
simulation directly to the data. 

About 8 X lO"^ upstream stop events were considered 
for this analysis, with 1.5 x 10^ events accepted after cuts. 
About ten times that amount of simulated decays were 
generated. 

As mentioned earlier, to first order energy loss is pro- 
portional to the amount of material a particle encoun- 
ters; in twist's planar geometry this was proportional 
to l/cos6'. Multiplying Ap by cos 9 removes this first- 
order angle dependence. Figure [7] shows the distributions 
of (Ap)(cos0), for both data and simulation; the simula- 
tion has been normalized to data by integrated counts. 
The most probable value was —28.4 ±0.1 keV/c for data 
and -29.65 ± 0.04 keV/c for MC; the FWHM width was 
155.9 ± 0.1 keV/c for data and 141.64 ± 0.04 keV/c for 
MC. The slight difference in most probable energy loss 
was compensated by the energy calibration. The differ- 
ence in the width of the distribution leads to a system- 
atic error in the measurement of the decay parameters of 
1.3xl0~^, as discussed in Sect. IVllEl and a correction 
was applied (Table UlI)) . Figure [5] shows the same distri- 
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FIG. 8: (Color online.) As Fig. [T] for large energy losses. 

Figure [9] shows the distributions of A9 for data and 
simulation; the simulation was normalized to the data as 
before. The most probable value was — 0.97± 0.02 mrad 
for data and -0.581 ± 0.007 mrad for MC; the FWHM 
width was 29.75 ± 0.02 mrad for data and 29.159 ± 
0.007 mrad for MC. The difference in the most probable 
A9 values between data and MC is very small compared 
to the approximately 10 mrad angular resolution and so 
did not affect the decay spectrum. Unlike the energy loss, 
the slight difference in the width of the A6 distributions 
had a negligible effect on the measurement of the decay 
parameters. Figure [TO] shows the same distributions but 
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FIG. 9: (Color online.) Scattering angle distributions for 
positrons across the target region (Al and PCs), showing 
events with small angle scatters and integrated over the fidu- 
cial region. 



FIG. 7: (Color online.) Energy loss distributions for positrons 
across the target region (Al and PCs), showing events with 
small losses and integrated over the fiducial region. 

butions but for events where (Ap)(cos6') < — 1 MeV/c, 
demonstrating very good agreement between data and 
simulation over several orders of magnitude. This was 
used to determine part of the systematic uncertainty due 
to the simulation of positron interactions (Sect. FVIIDIl . 
which was found to be less than 0.7x 10^* for both p and 



including events with large A9, again demonstrating very 
good agreement between data and simulation over many 
orders of magnitude. 

The upstream stops technique was also used to com- 
pare reconstruction efficiencies in data and simulation, 
as a function of momentum and angle. A simplified 
definition of reconstruction efficiency was used to avoid 
introducing artifacts at the edge of the fiducial region 
and other complications. Let rd{p, cos 9) be the num- 
ber of downstream tracks in a given (p, cos 9) bin for 
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FIG. 10: (Color online.) As Fig. [9] for large angle scatters. 



which an upstream track was also reconstructed, and 
let Td{p, cos 9) represent the total number of downstream 
tracks in that bin. Then the upstream efficiency was de- 
fined as eu{p,cos9) = ra/Td. Downstream efficiency was 
similarly defined. 

Figures [TT] and fT2l show the (Data-MC) differences in 
reconstruction efficiency, as functions of momentum and 
cos 9 respectively; the differences of upstream efficiencies 
are shown, and the downstream efficiencies were simi- 
lar. Each bin represents an average across the fidu- 
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FIG. 11: Difference between reconstruction efficiencies be- 
tween data and simulation (MC) within the fiducial region as 
a function of momentum, for the upstream half of the detec- 
tor. The difference for the downstream half is similar. 

cial region. The region around the upstream stops beam 
momentum of ~ 25 MeV/c has been excluded to avoid 
an artifact in the efficiency measurement caused by the 
beam positron phase space. The simulation and data 
show good agreement in reconstruction efficiency over 
the entire fiducial region, and the slight differences are 
independent of momentum and angle. Overall, we find 
that all but (5.7±0.2)xl0-'' [(5.12±0.05) x 10"'*] of the 
upstream stops events in the data [simulation] that con- 
tain a reconstructed downstream track within the fiducial 
region also contain a reconstructed track in the upstream 
half of the detector; similarly, all but (8.2 ± 0.2) x 10"'' 
[(7.87 ± 0.06) X 10"^] of the events that contain a recon- 




FIG. 12: Difference between reconstruction efficiencies be- 
tween data and simulation (MC) within the fiducial region as 
a function of cos 9, for the upstream half of the detector. The 
difference for the downstream half is similar. 



structed upstream track within the fiducial region also 
contain a reconstructed track in the downstream half of 
the detector. The differences in upstream and down- 
stream reconstruction efficiency give rise to a systematic 
error of less than 0.6x10"'', as discussed in Sect. IVIfBl 



VII. SYSTEMATIC ERRORS AND 
UNCERTAINTIES 

A. General Procedure 

A systematic uncertainty is an uncertainty in the decay 
parameters as the result of an unknown or unaccounted- 
for value or variation of some experimental parameter. 
A summary of the systematic uncertainties for this mea- 
surement is given in Table |TT] by category; these are de- 
scribed in more detail below. A systematic error repre- 

TABLE II: Summary of systematic uncertainties by category. 



Category 



Ap 



A5 



Chamber Response 
Energy Scale 
Positron Interactions 
Resolution 



0.00029 0.00052 
0.00029 0.00041 
0.00016 0.00009 
0.00002 0.00003 



Alignment and Lengths 0.00003 0.00003 

Beam Intensity 0.00001 0.00002 

Correlations with -q 0.00011 0.00001 

Theory 0.00003 0.00001 



Total 



0.00046 0.00067 



sents a known experimental bias rather than an uncer- 
tainty; where these errors were found to be significant, 
a correction was applied. Very small errors were simply 
included in the total systematic uncertainties. Table IIIII 
lists the corrections applied to the present measurement. 



10 



TABLE III: List of corrections to final decay parameter mea- 
surement, in units of 10"*. Uncertainties in the corrections 
are included in tfie systematic uncertainties in Table HIl 



Correction 



Ap 



AS 



Drift time maps -2.0 ± 2.9 +1.6 ± 5.2 

Momentum resolution +1.2 + 0.2 +1.3 ± 0.3 



Total 



-0.8 + 2.9 +2.9 + 5.2 



Several systematic uncertainties important to the mea- 
surement of PJ^£, were not studied for this measurement. 
The leading systematic uncertainty in the previous mea- 
surement of P^S, by TWIST [2J| arose from uncertainty 
in our knowledge of the muon beam emittance, and its 
impact on the muon depolarization as the beam entered 
the magnetic field. The depolarization of the muons due 
to the muon beam emittance was found to be 5-7x 10^'^ 
depending on the beam used, and the associated system- 
atic uncertainty was AP!^^ — 0.0034. The muons depo- 
larized further while at rest in the stopping target at a 
rate of (1.6 ± 0.3) x 10~^ f^s~^, leading to an associated 
systematic uncertainty of ^P^S, — 0.0012. We have no 
additional data to improve our knowledge of the muon 
beam used at the time. Furthermore, whereas this issue 
had a very large impact on the determination of PJ^S,, its 
impact on the measurement of p and 5 is negligible. 

Systematic uncertainties and errors were studied using 
the spectrum fitting technique described by Eqn. (fTO|) . 
A new simulated decay spectrum was produced by ex- 
aggerating an experimental parameter at the simulation 
or analysis stage, and this was fit against a simulation 
generated using the standard parameters. The fit results 
Ap and AJ, scaled by the amount of exaggeration, pro- 
vided a direct measurement of the systematic uncertainty 
associated with the experimental parameter. 

All systematic errors and uncertainties were deter- 
mined prior to revealing the hidden decay parameters 
assumed by the simulation. 



B. Chamber Response 

To determine the effect of inaccuracies in the 
GARFiELD-generated STRs, analyses using STRs derived 
directly from the data and the simulation were fit against 
the corresponding standard analyses. Data-derived and 
simulation-derived STRs for TWIST have only recently 
been developed. They were derived for data and simula- 
tion from the means of the time residuals (the difference 
between reconstructed hit time and measured hit time) 
as a function of hit position within the drift cell, during 
reconstruction using GARFiELD-generated STRs. These 
residuals were used to modify the GARFiELD-generated 
STRs, and the process was iterated until it converged. 
The data-derived STRs automatically account for varia- 
tions in temperature, foil positions, etc., within the de- 



tector, as well as compensating for some residual biases 
in the helix reconstruction. They lead to reduced x^ val- 
ues from the helix fitter. They also provide improved 
momentum resolution and a better match between data 
and simulation compared to that described in Sect. IVII 
Thus, we conclude that they are more correct for use in 
the analysis than the STRs generated by Garfield. 

The tests determined by how much the measured de- 
cay parameters were affected by the use of GARFiELD- 
generated STRs in the standard analysis. The difference 
between the effect in simulation and the effect in data 
demonstrated the amount by which the measured decay 
parameters in a standard data-MC fit would shift; a cor- 
rection was applied to the final measurement to account 
for this (Table [lTl|) . The uncertainty on this difference, at 
3xl0~^ for p and 5xl0~^ for J, represented the system- 
atic uncertainty due to inaccuracies in the STRs. This 
dominated the chamber response uncertainties. 

Other aspects of the chamber response studied in- 
cluded the uncertainties in the DC foil positions, the 
asymmetries in the decay reconstruction efficiency, and 
electronics-related wire time offsets. None of these rep- 
resented uncertainties greater than 0.6 x 10~* for either 
p or (5, and most were much smaller. 



C. Energy Scale 

To translate uncertainties in the energy calibration pa- 
rameters into uncertainties in the decay parameters, new 
spectra were created by applying energy calibration, with 
each calibration parameter exaggerated, to a standard 
simulation. These exaggerated spectra were fit against a 
spectrum created using the standard calibration to mea- 
sure the effects on the decay parameters, which were 
found to be 3xl0~* for p and 4xl0~* for 5. The sys- 
tematic uncertainty due to the energy calibration was the 
dominant uncertainty related to energy scale. 

Other aspects of the energy scale studied included the 
assumed behavior of the energy calibration with momen- 
tum, and errors in the shape of the magnetic field map 
used for analysis. The energy calibration was assumed 
to be independent of the positron momentum; using a 
calibration proportional to the momentum changes the 
decay parameters by less than 0.8 x 10^'' for both p and 
5. The simulated magnetic field used for analysis and 
simulation was compared against the measured magnetic 
field map; the differences were found to influence the de- 
cay parameters at the level of 0.7x10^* for both p and 
5. 



D. Positron Interactions 

The production of delta rays and bremsstrahlung were 
the most important discrete positron interactions for 
TWIST. The simulation of these was validated using the 
upstream muon stops described in Sect. IVTl above. The 
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rate of delta ray production in simulation was compared 
to that in data using the ratio Rs — N12/N11, where N12 
is the number of events with one upstream track and 
two downstream tracks, and A^n the number of events 
with one upstream track and one downstream track. We 
find Rs = (1.432 ± 0.003) x 10"^ ^^ ^he upstream stops 
data. Rs was also measured in several additional up- 
stream stops simulations with the delta ray production 
cross-section multiplied by various factors, as shown in 
Fig. [131 A linear fit to Rs as a function of the cross- 
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FIG. 13: Rs as a function of delta ray cross-section multiplier. 
Open circles are simulation; standard simulation corresponds 
to a multiplicative factor of 1. The linear fit was used to 
determine an effective multiplier for the data (dashed line). 
Error bars are shown but are vanishingly small. 

section multiplier was used to translate the Rs value 
from data into a difference in delta ray cross-sections be- 
tween data and simulation, conservatively assuming that 
all events of this topology are due to delta rays. This is 
not true in practice, since Rs does not go to zero in the 
simulation when the delta ray production cross-section in 
GEANT is set to zero (see Fig. [T3)) : however, the conser- 
vative assumption was sufficient for estimating the sys- 
tematic uncertainty. The effective delta ray production 
cross-section determined from this method was found to 
be approximately 18% lower in the simulation than was 
seen in data. 

The rate of bremsstrahlung production (Rb) in sim- 
ulation was compared to that in data by counting the 
number of through-going positrons whose change in mo- 
mentum was pds — Pus < ii^p) ^ 1 MeV/c), normal- 
ized to the total number of reconstructed events; here 
(Ap) represents the most probable value of Ap. We find 
Rb = (1.42±0.01)xl0"2^ As with the study of delta ray 
production, Rb was also measured in several additional 
upstream stops simulations with the bremsstrahlung pro- 
duction cross-section increased by various factors, as 
shown in Fig. [TD As before, a linear fit was used to 
translate the Rb value from data into a difference in 
bremsstrahlung cross-sections between data and simula- 
tion; the simulation was found to agree with the data to 
within 2%. 

The simulations produced with the delta ray or 
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FIG. 14: Rb as a function of bremsstrahlung cross-section 
multiplier. Open circles are simulation; standard simulation 
corresponds to a multiplicative factor of 1. The linear fit was 
used to determine an effective multiplier for the data (dashed 
line). Error bars are shown but are vanishingly small. 



bremsstrahlung cross-sections increased by a factor of 
three were fit against standard simulation. The sys- 
tematic uncertainty due to errors in the simulation of 
delta ray production was found to be 1.5x 10^"* for p and 
0.9x10""' for S; that due to errors in the simulation of 
bremsstrahlung production was less than 0.7xl0~^ for 
either parameter. 

Materials outside the sensitive region of the detector 
could scatter particles back inside, resulting in extra- 
neous hits that the track reconstruction must sort out. 
The number of backscattered positrons, normalized by 
the number of muons, was studied in the data and the 
simulation. A new simulation was produced with an ad- 
ditional plate of aluminum downstream of the detector 
as a backscattering source, and this was fit against the 
standard simulation. The effect on the decay parameters 
was less than 0.4x10^"*. 



E. Resolution 

Reconstruction resolution in momentum and angle 
smeared the decay spectrum. The RMS angle resolution 
was about 5-15 mrad, and the simulation agreed with 
the data to within 3 mrad. The RMS momentum res- 
olution was 0.040-0.120 MeV/c, and was smaller in the 
simulation by 0.007 MeV/c. Both types of reconstruc- 
tion resolution varied with energy and angle of the de- 
cay positron. Since the spectrum was inherently smooth, 
however, these resolutions did not significantly distort 
the decay parameters directly. The largest effect was at 
the endpoint, where the momentum resolution changes 
the shape of the edge; the 7 keV/c difference in resolution 
between simulation and data had a small effect on the 
energy calibration, which resulted in a distortion in the 
decay parameters. Resolution was measured in the data 
and MC as a function of momentum and angle using de- 
cays from upstream muon stops, by measuring the widths 
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of the energy loss and scattering distributions. New spec- 
tra were produced by smearing the reconstructed mo- 
mentum or angle of each event in the simulation before 
including it in the spectrum and calculating the energy 
calibration. The effect of the angle resolution on the de- 
cay parameters was found to be less than 0.4x10^"' for 
both p and 5. The effect of the momentum resolution on 
the decay parameters was 1.2x10^'* for p and 1.3xl0~^ 
for (5, and a correction was applied (Table Hill). 



F. Other Sources of Error 

A number of other possible sources of systematic error 
were studied and found to be small. 

The effect of the uncertainty in the chamber alignment 
on the decay parameters was tested by analyzing a sim- 
ulation using distorted chamber alignments, and fitting 
against the standard simulation. These were found to 
affect the decay parameters by less than 0.03x10""'. Er- 
rors in the transverse and longitudinal length scales of the 
detector translated directly into errors on the transverse 
and longitudinal components of the momentum; these 
were tested by distorting these momentum components 
for reconstructed decays before including them in a new 
spectrum. The effect was less than 0.3x10^*. 

New simulations were produced with exaggerated 
muon or positron beam rates, to test the sensitivity of 
having an incorrect pile-up rate in the simulation. This 
showed the systematic uncertainty due to the simulation 
of beam intensity to be less than 0.2x 10~^ for both beam 
components. 

As described above, the highest order of radiative cor- 
rections used for this analysis was 0{a^L). The 0{a^) 
radiative corrections [221 represented the theoretical er- 
ror in the 0{a'^L) corrections, and this determined the 
theoretical uncertainty for this analysis. A new simula- 
tion was produced with the 0{a^L) radiative corrections 
exaggerated, and this was fit to the standard simulation. 
The results show that the theoretical uncertainty in the 
measured decay parameters is less than 0.3x10"'*. 

During decay parameter fits, rj was normally fixed to 
the world average value of -0.0036 25]. The standard 
fit between data and MC was repeated with rj raised or 
lowered by one standard deviation (^77 = ±0.0069), giv- 
ing the correlations dp/drj = 0.0162, dS/dr] = 0.0015, 
and dPfj_£,/d7] = 0.0155. This led to an uncertainty in p 
due to the uncertainty in 77 of 0.00011. Future improve- 
ments in the knowledge of 77 can be used to reduce this 
systematic uncertainty directly. 

Since this measurement was finalized, development of 
the simulation and analysis software has continued for the 
analysis of additional data. During this process, an er- 
ror in the particle identification was found, which caused 
good events to be misclassified and discarded, affecting 
about 0.6% of events in the fiducial region. This distor- 
tion was primarily linear in |cos0|, which is a shape not 
refiected in any of the derivative spectra; furthermore, 



the distortion was the same in the analysis of both data 
and simulation, leaving the derived decay parameters un- 
affected. This software error has been fixed for future 
analyses, but has been listed here for completeness. 



VIII. RESULTS 

The values of p and 6 measured from fitting each data 
set with its corresponding simulation are listed in Ta- 
ble IIVI After taking weighted averages and applying the 
corrections in Table IIIIl we find 



p = 0.75014 ± 0.00017(stat) ± 0.00044(syst) ± 0.00011(7^), 

where the last uncertainty is due to the uncertainty in 77, 
and 

6 ^ 0.75067 ± 0.00030(stat) ± 0.00067(syst). 

Both results are consistent with the Standard Model val- 
ues of 3/4. These represent a factor of two improvement 
over the previous TWIST measurements [lO|, lll| . 



TABLE IV: Measured decay parameters from fits between 
data and simulation, for each data set. Uncertainties are sta- 
tistical. Each fit has 2463 degrees of freedom. 

Set p S x^ 

Mis-steered 0.75054 ± 0.00044 0.75066 ± 0.00077 2505 
PC5 stops 0.75112 ± 0.00079 0.74757 ± 0.00140 2434 
Stop i 0.74977 ± 0.00038 0.75081 ± 0.00067 2458 

Stop I A 0.75024 ± 0.00037 0.75073 ± 0.00066 2483 
High rate 0.75003 ± 0.00050 0.75047 ± 0.00088 2392 
Aperture 0.75019 ± 0.00045 0.75116 ± 0.00080 2555 
Stop I B 0.75042 ± 0.00049 0.74882 ± 0.00086 2468 



The typical correlation coefficient between p and S 
from the decay parameter fits is -f0.15 for the fiducial 
range adopted here. Correlations also exist between the 
systematic uncertainties in p and d. The contributions to 
a^g associated with the chamber response and positron 
interaction systematics are also positive and somewhat 
larger than the statistical contribution from the decay 
parameter fits. In contrast, the energy calibration sys- 
tematics for p and 6 are strongly anti-correlated, and 
provide the dominant contribution to a^g in the present 
measurement. Overall, we find that the correlation coef- 
ficient between the total uncertainties in p and 6 is —0.16. 

The decay parameter fits also determined a value for 
P^^, which can be converted to a value for P^^ by cor- 
recting for known depolarization effects. After apply- 
ing corrections corresponding to the p and 6 corrections 
in Table UTTl the fits give PJ ^ = 1.0025 ± 0.0004(stat). 
As explained in Sect. IVII A[ the systematic uncertain- 
ties important to the measure of P^^ were not revisited 
during this analysis, and this value should not be con- 
sidered a revised measurement. When differences in the 
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TABLE V: Results of a new global analysis of muon decay 
data, including the present measurements (90% C.L.). P^ — 
1 is assumed. Best fit values of selected decay parameters are 
also listed. 

Gagliardi ct al. [25] Present Analysis 
(xlO~^) (xlQ-^) 



TABLE VI: Limits on the weak coupling constants. (Limits 
on Ifffil and Ig^^l are from Ref. JH.) 



Qrr 


< 1.14 


< 0.96 


Qlr 


< 1.94 


< 1.38 


Qrl 


<44 


<42 


Qll 


>955 


> 955 


Blr 


< 1.27 


<0.64 


Brl 


< 10.9 


< 10.8 


a/A 


0.3 ±2.1 


0.1 ±1.6 


P/A 


2.0 ±3.1 


2.1 ±3.0 


a' /A 


-0.1 ±2.2 


-0.1 ±1.6 


P'/A 


-0.8 ±3.2 


-0.8 ±3.1 



P 


0.74959 ± 0.00063 


0.74964 ± 0.00035 


S 


0.74870 ±0.00114 


0.74997 ± 0.00065 


V 


-0.0036 ± 0.0069 


-0.0042 ± 0.0064 



two analyses are considered, the value of P^^ found here 
is consistent with the published TWIST measurement of 
P^^ = 1.0003 ± 0.0006(stat) ± 0.0038(syst) (H. 





Gagliardi et al. 


[25] Present Analysis 


\9rr\ 


< 0.067 


< 0.062 


\gU 


< 0.034 


< 0.031 


IqU 


< 0.088 


< 0.074 


\9Ir\ 


< 0.036 


< 0.025 


\glR\ 


< 0.025 


< 0.021 


\gL\ 


< 0.417 


< 0.412 


\g}iL\ 


< 0.104 


< 0.104 


\91l\ 


< 0.104 


< 0.103 


\gL\ 


< 0.550 


< 0.550 


IffLl 


> 0.960 


> 0.960 



The values of the Qei_i from this global analysis can 
be used in Eqn. ^ to place limits on the magnitudes 
of the weak coupling constants \g]n,U the exceptions are 
gYi^ and gf^^^ , which are determined more sensitively 
from inverse muon decay, e~v^ — > fJ-^Ve- The limits de- 
termined with this method are listed in Table IVll along 
with the values from the previous global analysis [25J . 
The present analysis represents a reduction of approxi- 
mately 16% in the limits for |<?f^| and Ig^^^l, and a 30% 
reduction for Ig)^^!. 



IX. GLOBAL ANALYSIS OF MUON DECAY 



X. CONCLUSION 



In 2005 Gagliardi et al. performed a global analysis 
of all available muon decay data 12511. in cluding earlier 
TWIST measurements of p and 6 |10l.lllll. That analysis 
has been repeated, incorporating the TWIST measure- 
ment of P^S, [24l] and the new p and S measurements 
presented here; the correlation factor of —0.16 between 
the TWIST p and 6 measurements has also been included 
in the calculation. All other input values are the same as 
in the analysis of Gagliardi et al. 

In brief, the global analysis used a Monte Carlo method 
similar to that of Burkard et al. [2^ to map out the 
joint probability distributions for 9 independent vari- 
ables, Qrr, Qlr, Qrl, Blr, Brl, a/ A, [i/A, a' /A, 
and P' /A. Each of these parameters is a bilinear com- 
bination of the weak coupling constants gj^. The decay 
parameters could then be written in terms of these in- 
dependent variables. Table |V] shows the results of this 
global analysis, as well as the results of the analysis of 
Gagliardi et al. 

The 90% confidence limits are given for the indepen- 
dent variables listed above, and global best-fit values of 
the decay parameters p, 5, and rj are given. The present 
analysis represents significant improvements in the lim- 
its on Qlr and Blr, and tightens several of the other 
limits. It is interesting to note that the global analysis 
significantly reduces the uncertainty in the value of p, 
from a total of 0.00063 to 0.00035. 



This new measurement of the muon decay spectrum 
is a factor of two more precise than previous measure- 
ments [i3, [m ■ It is consistent with Standard Model pre- 
dictions, placing more stringent limits on "new physics" 
in the weak interaction. New indirect limits on the value 
of P^^ can be obtained using these values of p and 8, 
in combination with the the measurement of Pn£,S/p > 
0.99682 at 90% confidence by Jodidio [13, Hf] and us- 
ing the constraint £.6/ p < 1 (required to obtain positive 
definite decay probabilities from Eqn. ^). Accounting 
for the correlation coefficient of —0.16 between the p and 
S uncertainties, we find 0.99524 < P^^ < £. < 1.00091 
at 90% confidence. This is a significant improvement 
over the previous indirect limit of 0.9960 < PJ^ < ^ < 
1.0040 [li|. 

The quantity Q^ = Qrr + Qlr represents the total 
probability for a right-handed muon to decay into any 
type of electron, a process forbidden under the Standard 
Model weak interaction. The new measurements of p 
and 6 lead to the new limits on Qrr, and Qlr shown 
in Table |Vl and hence to a new 90% confidence limit 
upper bound on the combined probability Q^ < 0.0024, 
a slight improvement over the limit of Q^ < 0.003 from 
the previous global analysis [25| . 

Under left-right symmetric models, p > 0.75 is forbid- 
den, so the general measurement of p can be converted 
into a 90% confidence limit lower bound within LRS mod- 
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els: p > 0.7493 (compared with p > 0.7487 from the p 
measurement previously published by TWIST). The re- 
lation p ~ |(1 — 2Cg) then gives a 90% confidence limit 
of ICgl < 0.022, a significant improvement over the limit 
of ICgl < 0.030 for the previously pubHshed TWIST value 
of p. 

The final phase of TWIST is in progress. Additional 
data are in hand, taken in 2006 and 2007, and new anal- 
ysis is underway. Further improvements to the simula- 
tion and analysis are being implemented, including the 
use of measured drift time tables and improved track re- 
construction algorithms, and the resulting better agree- 
ment in the momentum resolution between simulation 
and data. These improvements are expected to lead to 
additional factor-of-two reductions in the uncertainties 
on p and 6, providing another incremental improvement 



to searches for new physics. 

We would like to thank C.A. Ballard, S. Chan, 
J. Doornbos, B. Evans, M. Goyette, D. Maas, J. A. Mac- 
donald (deceased), T.A. Porcelli, N.L. Rodning (de- 
ceased), J. Schaapman, G. Stinson, V.D. Torokhov, 
M.A. Vasiliev, and the many undergraduate students 
who contributed to the construction and operation of 
TWIST. We also acknowledge many contributions by 
other professional and technical staff members from TRI- 
UMF and collaborating institutions. This work was sup- 
ported in part by the Natural Science and Engineering 
Research Council of Canada and the National Research 
Council of Canada, the Russian Ministry of Science, and 
the U.S. Department of Energy. Computing resources for 
the analysis were provided by WestGrid. 



[1] W. Fetscher, H.-J. Gerber, and K. Johnson, Phys. Lett. [17] 

B 173, 102 (1986). 
[2] W. Fetscher and H.-J. Gerber, in ^. 

[3] L. Michel, Proc. Phys. Soc. A 63, 514, 1371 (1950). [18] 

[4] C. Bouchiat and L. Michel, Phys. Rev. 106, 170 (1957). 
[5] T. Kinoshita and A. Sirlin, Phys. Rev. 108, 844 (1957). 
[6] W.-M. Yao et al., J. Phys. G 33 (2006). [19] 

[7] A. Davydychev, K. Schilcher, and H. Spiesberger, Eur. [20] 

Phys. J. C 19, 99 (2001). 
[8] P. Herczeg, Phys. Rev. D 34, 3449 (1986). [21] 

[9] S. Profumo, M. Ramsey-Musolf, and S. Tuhn, Phys. Rev. 

D 75, 075017 (2007). [22] 

[10] J. Musser et al. (TWIST Collaboration), Phys. Rev. Lett. 

94, 101805 (2005). [23] 

[11] A. Gaponenko et al. (TWIST Collaboration), Phys. Rev. [24] 

D 71, 071101(R) (2005). 
[12] C. Oram, J. Warren, G. Marshall, and J. Doornbos, Nucl. [25] 

Instrum. Methods 179, 99 (1981). 
[13] R. Henderson et al., Nucl. Instrum. Methods A 548, 306 [26] 
(2005). [27] 

[14] J. Hu et al., Nucl. Instrum. Methods A 566, 563 (2006). [28] 

[15] Y. Davydov et al., Nucl. Instrum. Methods A 461, 68 

(2001). 
[16] Vector Fields Ltd., Computer code OPERA-Sd (2002). 



The TWIST Data Acquisition System at TRIUMF (Mon- 
treal, QC, 2003), TRIUMF preprint TRLPP-04-16 (un- 
published) . 

R. Veenhof, GARFIELD: Sim,ulation of gaseous detec- 
tors, CERN (2003), version 7.10, CERN program Library 
long writeup W5050 (unpublished). 
A. Arbuzov, Phys. Lett. B 524, 99 (2002). 
A. Arbuzov, A. Czarnecki, and A. Gaponenko, Phys. 
Rev. D 65, 113006 (2002). 

A. Arbuzov and K. Melnikov, Phys. Rev. D 66, 093003 
(2002). 

C. Anastasiou, K. Melnikov, and F. Petriello, J. High 

Energy Phys. 09, 014 (2007). 

G. Lutz, Nucl. Instrum. Methods A 273, 349 (1988). 

B. Jamieson et al. (TWIST Collaboration), Phys. Rev. 
D 74, 072007 (2006). 

C. Gagliardi, R. Tribble, and N. Williams, Phys. Rev. D 
72, 073002 (2005). 

H. Burkard et al.. Physics Letters B 160, 343 (1985). 
A. Jodidio et al., Phys. Rev. D 34, 1967 (1986). 
A. Jodidio et al., Phys. Rev. D 37, 237 (1988), erratum 
to fH. 



